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Abstract 

We study the virial coefficients Bk of hard spheres in D dimensions by means of Monte-Carlo in- 
tegration. We find that B5 is positive in all dimensions but that B§ is negative for all D > 6. For 
7 < k < 17 we compute sets of Ree-Hoover diagrams and find that either for large D or large k the 
dominant diagrams are "loose packed" . We use these results to study the radius of convergence and the 
validity of the many approximations used for the equations of state for hard spheres. 
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1 Introduction 

The question of the possible negativity of virial coefficients in the low density expansion 

^=„ + i>/ a) 

a k=2 

of the system of hard spheres with diameter a in D dimensions specified by the two body pair potential 

™-{T tit: < 2 > 

has been an unresolved problem of outstanding importance since it was first proposed by Temperley in 
1957. For dimensions D < 5 all currently available information is summarized in Table ^ where it is seen 
that all E>k are positive. However it was first shown [2] in 1964 for D > 8 that B4 is indeed negative. The 
best available current results for B4 as a function of dimension are shown in Table [21 

Most of our intuition and physical insight into the low density (fluid) phase of hard spheres in 3 dimensions 
comes from the 8 term virial expansion of Tabled Over the years this data has been used to produce many 
approximate equations of state 011311112312112212312112312111231211123123122 These 
approximates all incorporate the feature of positive virial coefficients and they all have the feature that they 
have a radius of convergence which is greater than the packing fraction r\j — 0.49 at which freezing has been 
seen to occur in computer experiments jMJ E3 El where the packing fraction 77 is related to the density p by 
V = B2p/2 D ~ 1 ■ This analyticity at the freezing density is incorporated into most phenomenological theories 
of freezing |3oi 1371 I38| as a homogeneity or mean field approximation which ignores the fluctuations at phase 
coexistence between the fluid and solid phases [3"§| . 

However, if there are negative virial coefficients for hard spheres in D = 3 then no conclusion on the 
radius of convergence of the virial expansion based on Tabled can be considered as reliable. 
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e-mail: Nathan.Clisby@stonybrook.edu and mccoy@insti.physics.sunysb.edu 
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Tabic 1: Numerical values of the virial coefficients B^jB\ 1 for k = 3, . . . , 8 for D = 2, 3,4, 5. 
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Table 2: Exact and numerical results for B 2 , B3, and B4 for 2 < D < 12. 
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The most striking effect of negative viral coefficients will occur if the signs oscillate with some period 
as k — > oo because this will give a radius of convergence which is not on the positive real axis. If this 
radius is less than the freezing density then it will be impossible to reliably learn anything about the freezing 
transition from a knowledge of a finite number of virial coefficients. 

It is thus most significant that the sixth and seventh virial coefficients for parallel hard cubes 0H] were 
shown to be negative in 1962, and even more important that for an exactly solved hard squares model 41 
and the hard hexagon model [421 1431 14*4^ the radius of convergence is limited by a singularity at complex 
density thus resulting in virial coefficients that oscillate in sign. 

The purpose of this paper is to extend the results of Tables ^and E] and examine as closely as possible the 
question of whether or not the virial coefficients B^ of the hard sphere gas in D dimensions have negative 
virial coefficients. The method we shall use is Monte-Carlo evaluations of the integrals in the Ree-Hoover 
expansion. In Section |2] we review the formalism of the Ree-Hoover expansion to establish our notation. In 
Section |3] we compute the virial coefficients B$ and Bq for dimensions up to D = 50. We find that B§ is not 



2 



monotonic but is in fact always positive. More importantly we find that B$ is negative for all D > 6. Our 
numerical results are given in Table |3J For higher virial coefficients the number of contributing diagrams 
rapidly increases. Consequently in this study we restrict our attention to various classes of diagrams which 
are studied in Section 0] where we are able to determine the class of diagrams which are dominant for large 
k. In Section [S] we use our results to form estimates of the radius of convergence to the virial series and we 
conclude in Section with an evaluation of the various approximate equations of state for hard spheres. 

Table 3: Numerical results for B^/B^, B5/B2, and B§jB\. The underline indicates the position of the local 
minima and maxima. Values for each coefficient for D = 3, and for Bq in D = 4, 5 are taken from Tabled 
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2 Ree— Hoover expansions 

The original graphical expansion for the virial series is due to Mayer and Mayer |45j . in which each bond 
represents the function 

/(r)=exp(-C/(r)/fc B T)-l (3) 

where r is the distance between the two vertices. A useful re-summation was performed by Ree and Hoover [§J 
I46| by introducing the function 

/» = 1 + /(r) = cxp {-U(v)/k B T) (4) 

and then expanding each Mayer graph by substituting 1 = / — / for pairs of vertices not connected by / 
bonds. This method was previously used by Percus and Yevick |47| in comparing the exact values of the 
fourth and fifth virial coefficients with coefficients obtained from the Percus- Yevick equation, and by Percus 
0Hj in discussing the derivation of the Percus- Yevick equation. 
The fourth virial coefficient may then be written as 

*=-\w-\\z[-\n=\*-\\ He-Id (5) 



3 



where the first expression is the expansion in Mayer graphs, the second is the expansion in Ree-Hoover 
graphs with the / bonds indicated by dotted lines and the third shows the equivalent Ree-Hoover graphs 
with the / bonds indicated by solid lines. In the second expression the graph with no / bonds is represented 
by 0. In the case of hard spheres, the potential is given by Eq. [21 so /(r) and /(r) are particularly simple: 



The virial coefficient B k is given in terms of the Ree-Hoover diagrams of k points of which m are the end 
points of / bonds S k [m, i] and a combinatorial factor C k [m, i] as 

B k =J2 B k[m,i] (8) 

with 

B k [m,i] = C k [m,i]S k [m,i] (9) 
where the index i labels the graphs in the class with fixed k, m. The combinatorial factor is expressed as 

C k [m,i] = --TT -s k [m,i]a k [m,i] (10) 

with 

s k [rn,i] = k\/#AutSk[m,i] (11) 

where §AutS k [m, i] is the cardinality of the automorphism group of the diagram S k [m, i] and a m ,/[fc] is the 
"star content" as defined by Ree and Hoover g^I ■ We will let S k [m, i] represent both the Ree-Hoover graph 
and the value of the corresponding integral. The k dependence of C k [to, i] is calculated by using the following 
relation for the star content |46| . 

~a k [m,i] = (-l)*" 1 ^ -2)2 fc _i [m,i] (12) 

The diagram S k [m,i] has k — to points that are connected to all other points by / bonds and are therefore 
indistinguishable, leading to the relation 

j AutSfc [m, i] = (k — iri)l$AutS m [in, i] (13) 

to obtain for k > m 

C k [m,i] = (-l) k ^-^ 2 ( k ~\)c m [m,i] (14) 
\m — 1/ 

where we note that for the complete star diagram B k [Q, 1] we have C k [0, 1] = 1/fc. 

For arbitrary D the number of Mayer graphs grows asymptotically as k — > oo as |49| 

ofc(fc-l)/2 

N(k) ^ (15) 

However, when the Ree-Hoover transformation is made many diagrams have zero star content and hence do 
not contribute to the virial coefficient. From Table 21 it appears that the ratio of contributing Ree-Hoover 
diagrams to the number of Mayer diagrams is bounded below, and so it is reasonable to suppose that 

lim N RH (k)/N(k) > (16) 

k — >OG 



In addition Ree-Hoover graphs may be zero for geometrical reasons. The number of non-zero graphs 
for D = 2 are taken from HOEli an d listed in Table where the value in parentheses excludes diagrams 
found to be negligible but which were not proven to be zero. For D = 1 only one graph in the Ree-Hoover 
expansion is non-zero for each k, namely the complete star, but for D > 2 it is an open question as to how 
many non-zero graphs there are at order k. 
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Table 4: Number of contributing Ree-Hoover and Mayer diagrams as a function of order. The number of 
Ree- Hoover diagrams with non-zero star content for k — 9 is a new result, while the other values are taken 
from 01111 Ej. 
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1 


1 
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0.500 


0.410 


0.365 


0.366 
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2.1 Close packed diagrams 

The / graph form used by Ree and Hoover (e.g. in 46 ) has the property that a single graph represents 
the sequence of diagrams B^[m, i] with m fixed where the (k + l)th order diagram is obtained from the kth 
by adding an additional point connected to all other points by an / bond. For example, all complete star 
diagrams of order k may be represented by the same symbol 0. We call such graphs with m < k close packed 
because no two points can be further apart than 2a. 

We find from Eq. ^]that the sign of Bk[m, i] for fixed m, i is independent of k. We further obtain from 
Eq. ^]for k large and m, i fixed that 

C k [m,i] = C m [m,i] ( +0(k m - 2 )) (17) 

2.2 Loose packed diagrams 

For diagrams in the class Bk[k, i] there are no points which are connected to all other points by / bonds and 
therefore there exist sequences of diagrams which have the property that the size of the configuration grows 
as k — > oo. We refer to this class of diagrams as "loose packed". For these graphs the / bond notation is 
more convenient. 

The simplest loose packed diagram in Bk[k,i] is the simple ring of k of the / bonds. We denote this 
diagram by the symbol R. More generally we consider graphs where a point is replaced by a diagram. We 
call such diagrams insertion diagrams and when we wish to make the type of insertion visually apparent we 
use the notation R(-) to represent a ring with a point replaced by a diagram. We also find it useful to label 
insertion diagrams by R n! z[fc] where the index n is the number of points in the inserted diagram of k total 
points and I labels the diagrams with given k and n. All four and five point insertions were found by starting 
with all four and five point Mayer graphs, and adding an extra point that is connected to two points by / 
bonds. A canonical labeling for this graph was found using "nauty", a program due to B. McKay |5()j . and 
matched with the five and six point graphs of Ree and Hoover jS] to find the star content. All insertions 
with non-zero star content are given in Table [S] along with the size of their automorphism group and over 
all combinatorial factor. From Ea. ll4l we find the insertion diagrams have the alternating sign property that 
(—l) k R n! i[k] has a sign which depends only on n and I but is independent of k. 

Multiple insertions are also possible when k is sufficiently large. We denote such a graph with n insertions 
as R{™^ j.j {„ i } where the subscripts indicate the types of insertions. For Bq we see in Table [S] that 
there are three graphs which may be interpreted as being composed of two 4-point insertions. 

We give in Tables EHH1 the Ree-Hoover graphs and their associated combinatorial factors for k = 4,5,6. 
We here give both the representation of the graphs in terms of / and / bonds. For the class of loose packed 
graphs in B^ [fc, i] we indicate the interpretation in terms of insertion diagrams with either multiple insertions 
or with a ring of one point and two bonds. In three cases the identification is not unique. 

The labeling index i in Bk[m, i] is chosen such that in D = 2 the magnitude of the contribution decreases 
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Tabic 5: Four and Five point insertions 



Diagram 


Label 


Group size 


Lowest order 
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R5,5 


2 


6 


3 


-3(fc- l)/2 
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with increasing i. If the diagram is identically zero for D = 2 then the ordering obtained from I? = 3 is used 
when possible. 



Table 6: Ree- Hoover diagrams for B4. 



Label 


s fe [m, i] 


a k [m,i] 


C'k[m,i] 


/ form 


/ form 


Insertion 


£4 [0,1] 


1 


-2 


2/8 









B 4 [4, 1] 


3 


1 


-3/8 








□ 


R[4] 


Table 7: Ree- Hoover diagrams for £? 5 . 


Label 


Sfc [m,i] 


a k [m,i\ 


C k [m,i] 


/ form 


/ form 


Insertion 


£5 [0,1] 


1 


-6 


6/30 









£5 [4, 1] 


15 


3 


-45/30 












£5 [5,1] 


30 


-2 


60/30 










R 4 , 2 [5] 


£5 [5, 2] 


12 


1 


-12/30 










R[5] 


S 5 [5,3] 


10 


1 


-10/30 
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R4,l[5] 
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Table 8: Ree-Hoover diagrams for B 6 . For diagrams B 6 [6, 14], Bq[6, 16], and £? 6 [6, 17] the assignment of 
insertion diagram labels is not unique and both possible assignments are shown. 



Label 


s k [m,i] 


a k [m,i] 


C k [m,i] 


/ form 


/ form 


Insertion 


Be [0,1] 


1 


24 


-24/144 









Be [4, 1] 


45 


-12 


540/144 








Be [5,1] 
Be [5, 2] 
Be [5, 3] 


180 
72 
60 


8 
-4 
-4 


-1440/144 
288/144 
240/144 








B 6 [6,1] 
B 6 [6, 2] 
B 6 [6, 3] 
B 6 [6, 4] 
Be [6, 5] 
B 6 [6, 6] 
B 6 [6, 7] 
B 6 [6, 8] 
B 6 [6, 9] 
B 6 [6, 10] 
B 6 [6, 11] 
Be [6, 12] 
Be [6, 13] 
B 6 [6, 14] 

B 6 [6, 15] 
B 6 [6, 16] 

B 6 [6, 17] 

Be [6, 18] 


360 
180 
60 
60 
180 
90 
45 
360 
360 
60 
15 
180 
360 
90 

90 
180 

15 

10 


3 
-2 
1 

-6 
-5 
-4 

4 
-1 
-2 

4 
16 

3 

1 

-2 

-1 
1 

1 

4 


-1080/144 
360/144 
-60/144 
360/144 
900/144 
360/144 
-180/144 
360/144 
720/144 
-240/144 
— 240/144 
-540/144 
360/144 
180/144 

90/144 
-180/144 

-15/144 

-40/144 




o 

W 

Q 

e 
e 

<x> 
o 


R 6 , 5 [6] 
R 4 , 2 [6] 

R[6] 
R 5 ,9[6] 

R{4,2}{4,2} [^1 
R 6 ,4[6] 

R5,s[6] 
Rs,7[6] 

R{4,1}{4,2} [6] 

or R5,6[6] 
Rs,3[6] 

R4,l[6] 

or R B ,i[6] 

R {4,1}{4 : 1}[ 6 ] 

or R 5 , 2 [6] 
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3 Monte-Carlo Calculation of Bk for D > k — 1. 

When D > k — 1 the integral Sk[m, i] can be written in a form where the dimension D appears as a simple 
power in the integrand. Hence the Monte-Carlo procedure can simultaneously calculate the given Ree-Hoover 
diagram in an arbitrary set of dimensions, including non- integer values. The key advantage of this method 
over that of [5j used in Section 0] is that we obtain fast convergence for high dimensions. We first present 
the details of the Monte-Carlo method and then the results for B4, B$, and B§. 

We believe that all Ree-Hoover diagrams are non-zero for D > k — 1 as it is possible to obtain a 
configuration for any diagram. We can see this by starting with a configuration where the distance between 
each point is exactly one, and then since each bond can be independently varied in length by a small amount 
we are able to satisfy the constraints imposed by / and / bonds of any graph. 



3.1 Monte-Carlo method 



In order to do the Monte-Carlo integration we need an appropriate measure for calculating Ree-Hoover 
diagrams in D-dimensional Euclidean space, where we wish to integrate out coordinates to leave ourselves 
with the lowest dimensional integral possible. The integrand is a product of / and / functions; for central 
potentials this leaves only the inter-particle distances as appropriate degrees of freedom. These form an 
independent set of k(k — l)/2 coordinates provided that D > k — 1, where k is the order of the diagram 
and hence the number of points in the configuration. After taking the infinite volume limit, for an arbitrary 
diagram Sk[m, i] we thus need to calculate integrals of the form 



/= []jd D r i F(r ij )=[l[n D -A [Y[da tJ [V({ aiJ })} D - k F( 
J i=l \i=l ) J i<j 



(18) 



27r D / 2 /r(,D/2), <n 



dji, V is the volume of the parallel-piped defined by these 



where flu—i 

distances in K , and F(rjj) is an arbitrary function of the inter-particle distances. As shown for example 
in |51| V may be expressed by the Cayley-Menger determinant: 



v({ ai }) = 



(~l) fc 

ofc-1 






1 


1 


1 


1 


1 





a 12 


013 ■ 


■ 0,1k 


1 


a2i 





a23 • 


■ a-ik 


1 


031 


032 


• 


■ a-ik 


1 


Ofcl 


a/c2 


Ofc3 ' 


■ 



(19) 



An overall scale factor can be taken out of Eq. ^| by taking 012 to be the largest value, and then setting 
a ij = a «j/ a i2j so that 



n f d < 

wjo JO 



r)n Jk-l)D/2-l 



^El.fl'.})^) 



(20) 



A formula similar to Ea. 1181 was previously obtained by Percus [52] . in which the integration of p vectors in 
dimension D is reduced to an integration of p vectors in p dimensions. 

For the hard sphere fluid, the integrand is either zero or ±1. In order to calculate Sk[m,i] with m f 
bonds and [k(k — l)/2 — m] f bonds, we may then proceed as follows. Generate a set of [k(k — l)/2 — 1] 
values uniformly distributed between and 1. Partition these values in two sets, where the largest m — 1 
values along with 1 represent / bonds and the other [fc(fc — l)/2 — m — X] values are / bonds, and then 



8 



randomly assign these values to edges in the diagram. Check that these values of the edge lengths squared 
can be embedded as a simplex in and if this is the case one can calculate the volume of the simplex 

and perform the an integral, making a contribution to the Monte-Carlo integral of 

i (* m {k - 1)D/2 -a m ^ 1)D/2 ) (ffn*-*) V ^ s PI) 

where a m and a m +i are the mth and (m + l)th largest values of respectively. 

We have used this procedure to compute B4, B 5 , and Bq. For B4 and B 5 500 batches of 5 x 10 6 con- 
figurations were used, while for Bq 20 batches of 2 x 10 9 configurations were generated. Uncertainties were 
calculated using 

(22) 

where there are q independent batches with value Ij. 



Error = 



J2 (< h > 

3=1 



<i>f 



<?(<?- 1) 



3.2 Results 

The values of the individual contributions Bk[m,i] to the virial coefficients Bk for k = 4, 5,6 are given in 
Appendix El in Tables EHH for many values of D < 50. The values of the virial coefficients B^/B^ 1 for 
D < 50 are given in Table 

The most important feature of the virial coefficients in Table |3 are the sign changes in B4 and Bq. A 
previous estimate for the dimension at which B4 becomes negative of D w 7.8 was obtained by Luban and 
Baram ^H] by means of a linear interpolation of the results of Ree and Hoover [2] , and also by an independent 
calculation. We further note the existence of a local minimum in B§ which confirms the tentative prediction 
of Loeser et al. jS3J. To obtain accurate values for these zeros, minima and maxima as a function of D 
the data were fitted with cubic splines, and all of these values are summarized in Table EB Note that the 
error estimates were made by fitting, for example, B4 + AB4 and B4 — AB4 with cubic splines to obtain 
a confidence interval for where B4 becomes negative. A least squares fit was not used because the error 
in values at different dimensions are not independent, as they are calculated simultaneously and not as 
independent samples. 

It is most important that Bq becomes negative at a lower dimension than B4. If this trend continues, 
we may expect that for £>2fc the dimension at which the coefficient becomes negative will decrease for higher 
orders and it is not unreasonable to expect that B% will be negative for D — 5 and possibly even D = 4. In 
addition, even though B$ is always positive, the existence of a local minimum at D = 6.87 shows that the 
increasing number of contributing diagrams results in more complex dimensional dependence. 



Table 9: Maxima and minima for B^/B^, B^jB\, and B§jB\. 



B4/BI 


B5/B2 


B 6 /Bl 


Becomes negative 
D\ e9 = 7.7320(4) 
Local minimum 
D = 10.7583(2) 


Local minimum 
D = 6.87 

Local maximum 
D = 8.31 


Becomes negative 
DJ 69 = 5.30(2) 
Local minimum 
D = 8.942(2) 



To obtain further insight into the structure of the virial coefficients we examine the individual contribu- 
tions in Tables IT21 - I14I There we see by examining the dependence of the contributions on D that by the 
time D = 50 the Ree-Hoover ring diagram R is several orders of magnitude greater than all other diagrams. 
We therefore make the following conjecture: 
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Conjecture 1 



lim B k (D)/R{D) = 1 

D— >oo 



(23) 



In particular it follows from this conjecture that for each k 

(-l) fe " 1 S fc (L>) > for sufficiently large D. (24) 

A similar conjecture in terms of Mayer diagrams was made by Frisch and Percus |54j in the course of their 
examination of Mayer diagrams at high dimensionality. 

We find further for D = 50 that not only does the Ree-Hoover ring diagram dominate but that the largest 
diagram in each B k [m, i] has the property that if m > m! then 

Bk[m, i] max >> B k [m',i] max . (25) 

We refer to this ordering as the principle of "loose packed dominance" . 

We also note that within the loose packed class B k [k, i] that the ordering of the magnitude of the diagrams 
at D — 2 and at D = 50 is drastically different. For example we note that the diagram R41 which vanishes 
identically in D = 2 is the second largest diagram for D = 50. More generally we find for D = 50 that the 
contributions in Bq[6, i] are ordered in magnitude as 

R > R 4 ,l > R4,2 > R 5,4 > R5,5 > R5,3 > R5,2 = R { 4 ,1},{4 : 1} 

(26) 

where we use the labeling of insertion diagrams in Tables [S] and |H1 We make the observation that these 
insertion diagrams are the largest in the class B k [k, i], and refer to this as "insertion graph dominance". 

4 Monte-Carlo Calculation of High Order Diagrams 

The observation and conjecture of the dominance of the Ree-Hoover ring diagram at large dimensions demon- 
strates that negative virial coefficients in the hard sphere fluids is a common occurrence. However it says 
nothing about the minimum dimension at which the coefficient B k will become negative for fixed k nor 
does it give any indication of whether or not virial coefficients of odd order can ever be negative. Ideally 
this question should be studied by computing further viral coefficients beyond those in Tabled for large k. 
Here we will begin this program by examining selected classes of diagrams B k [m, i] for values of k up to 17. 
In particular we examine the conjecture that the features of loose packed dominance and insertion graph 
dominance found for large D and fixed k holds also for fixed D and large k. 

4.1 Method 

The Monte-Carlo procedure used for this study is the same as that used by Ree and Hoover (e.g. in 9 ). 
Random configurations are generated by placing one point at the origin of some arbitrary coordinate system, 
then the next is randomly placed within a Z?-dimensional unit sphere centered on the first point, and so 
on. The final configuration is then tested to see if it satisfies the requirements imposed on distances by the 
diagram being calculated. For each diagram we used batches of 10 6 configurations, with enough batches to 
give accuracy better than two percent up to a maximum of 10000 batches, and errors were calculated using 

4.2 Results 

The evaluations of the following diagrams are presented in Tables I15H22I of Appendix [5] f° r values of k 
up to 17 and values of D up to 7; i?fe[0, 1], Sfc[4, 1], B fc [5, 1], R[fc], R4. 2 [fc], R4 j i[fc], R 5j5 [fc], and the "pinwheel 
diagram" in B k [k — 1, 1] obtained by adding one point which is connected by / bonds to all points in the 
Ree-Hoover ring R[fc — 1]. We note that for k — 5 the pinwheel is .85 [4, 1] and for k — 6 it is B e [5, 2]. 
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To investigate the conjecture of loose packed dominance for fixed D and large k we use the tables of 
Appendix 151 to find for each k and D which diagram of this set is the largest. This is shown in Table HOI 
We also use Tables ITBI and IT51 to compute the ratio R/B).[Q, 1]., and display the results in Table ITTI From 
this table we see that for D > 3 the Ree-Hoover ring R quickly dominates the complete star diagram as k 
increases. However, for D = 2, even though the ratios are monotonically increasing for k > 6 the ring has 
not dominated the complete star at k — 17. If we extrapolate the increasing ratios in Table ITTI for D = 2 we 
can estimate that the ring will be larger for approximately k ~ 22. 



Table 10: Largest diagrams 





Dimension 


k 


1 


2 


3 


4 


5 


6 


oo 


4 




















R 


5 

















R 


R 


6 














R 


R 


R 


7 











R 


R 


R 


R 


8 








RK-) 


R 


R 


R 


R 


9 








R«W 


R(<K) 


R 


R 


R 


10 








R(O0 


RM» 


R 




R 


11 








R(4>) 


R M» 


R 




R 


12 








R(O) 


R«W 


R«W 




R 


13 








R«W 


R(^) 


RW» 




R 


14 







R(<K> 


R«W 


RW» 




R 


15 







R«W 


R«N 


R(-O0 




R 



Table 11: R/0 



fc 


D = 2 


D = 3 


15 = 4 


D = 5 


15 = 6 


4 


-0.0299(2) 


-0.0942(4) 


-0.1964(8) 


-0.3405(7) 


-0.535(2) 


5 


0.0238(2) 


0.1143(7) 


0.362(3) 


0.911(6) 


1.97(1) 


6 


-0.0230(4) 


-0.184(2) 


-0.91(1) 


-3.47(5) 


-11.7(2) 


7 


0.0237(5) 


0.337(7) 


2.92(6) 


18.8(5) 


85(2) 


8 


-0.0277(6) 


-0.73(2) 


-10.4(3) 


-100(3) 


-780(30) 


9 


0.0350(8) 


1.67(5) 


44(1) 


710(30) 


11000(2000) 


10 


-0.043(1) 


-4.1(1) 


-177(7) 


-5100(800) 




11 


0.056(2) 


10.3(3) 


860(80) 






12 


-0.076(2) 


-26.0(7) 


-3000(800) 






13 


0.102(3) 


60(5) 








14 


-0.133(5) 


-160(30) 








15 


0.18(1) 










16 


-0.24(2) 










17 


0.39(5) 











Further information is obtained from the ratios of -Bfc[4, 1]/Bfc[0, 1] and Bk[5, 1]/Bk[0, 1] plotted in Figs. ^ 
and[5]of Appendix IO In these figures it is seen that the ratios -Bfc[4, l]/Bk[0, 1] and Bk[5, l]/£?fc[0, 1] increase 
linearly with k. Thus while for the values of k studied the ratios only rarely exceed one it appears that if we 
extrapolate to sufficiently large k that the inequalities in Eq. [23 will be satisfied for rn,m' — 0,4,5. 
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The relative size of diagrams in the class Bk[k, i] and Bk[k — 1, i] is studied by comparing the Rce-Hoover 
ring diagram R[k] of Table IT%1 with the corresponding pinwhccl diagram of Table l2~21 It is clear that in each 
dimension that for sufficiently large k the pinwheel diagram will vanish, and although this can't be shown 
numerically one can see the ratio of the pinwheel to the ring R[fc] decreases rapidly as k increases. 

We conclude this discussion of results by examining the relative magnitudes of insertion diagrams in 
Bk[k,i]. In Figs. andElof Appendix IH1 we plot the ratios R4,2/R, R4,i/R, and R5,5/R as a function of k 
for various D. We see that for sufficiently large k all three ratios have a linear increase with respect to k. We 
also note that even though the ratio R4.1/R4.2 is small in low dimensions, it approaches a non-zero constant 
as k — > 00, and for sufficiently large D the ratio is in fact greater that unity. 

The examination of all available data supports the following conjectures 

Conjecture 2 The loose packed diagrams in Bk[k, i] dominate all other diagrams for large k when D > 3. 

In connection with this conjecture we note that the three largest diagrams for Bq in D = 5,6 are Rs,5, R44, 
and R and that the sum of these three largest diagrams shares the property with Bq itself of changing sign 
between D — 5 and D = 6. We restrict the conjecture to dimensions D > 3 because at this stage there is 
not yet enough data to support the case for D = 2. 

Conjecture 3 The insertion diagrams all have the same exponential rate of growth as the Ree-Hoover ring. 



5 Estimates of the radius of convergence 

The dominance of loose packed diagrams for fixed D > 3 and large k may be used to discuss the question 
of sign change in the virial coefficients and the radius of convergence of the virial expansion. Of greatest 
importance is the relation of the radius of convergence to the packing fractions 77/ at which freezing occurs 
which have been numerically determined as 77/ = 0.49 in D = 3 [35II2II, Vf = 0-31 in D = 4 and r/f = 0.19 
in D = 5 as obtained from |55l 156) . 

It may be that there are sequences of loose packed diagrams that grow faster than the ring, but for low 
order the largest diagrams can be characterized as insertion diagrams and these appear to have the same 
exponential rate of growth as the ring. Hence we will concentrate on the Ree-Hoover ring R. We first note 
that the absolute value of this diagram must be strictly less than the absolute value of the corresponding 
Mayer diagram where all / bonds are replaced by unity. The resulting Mayer ring diagram has been long 
ago evaluated (S7| in terms of Bessel functions of the first kind J D / 2 (x) and thus we obtain 



(-l^Rm/B. 



k-l 



< 



(fc-l)(27T) 



kD/2 



2kn 



k-2 
D-l 



dx x 



D-l 



Jd/2(x) 



,D/2 



(27) 



The large k behavior of this integral is easily obtained by steepest descents by noting that the maximum 
value of Jd/2{x)/x d / 2 occurs at x = 0. Thus we have as k — ► 00 



(-l^RM/B*- 1 < 



(fc-l)(l + £>/2) g / 2 
k 1 + D / 2 Y{\ + D/2) 



,fc-2 



(28) 



and hence R[k] satisfies the bound 

|R[A + l]/(B 2 R[fc])| <2. (29) 

This leads to a packing fraction at the radius of convergence of the sum of Ree-Hoover diagrams of r\ r h = 2~ D 
which is substantially greater than the lower bound 0.145/2 15 of Lebowitz and Penrose but which is still 
much smaller than the freezing densities r\ t. This value of r]^ is of the same order of magnitude as the lower 
bounds on the termination density due to Baram and Fixman which are r\ t > 0.25838 for D — 2 and 
rjt > 0.12681 for D = 3. This is compatible with the pressure being analytic for positive values of 77 less than 
rjf because the alternations of sign of R[fc] puts the leading singularity on the negative r\ axis. 
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However the ratios \R[k + 1] / (i?2 R [fc] ) | obtained from Table H3I are substantially below the bound Eg. 1291 
In fact we see from each D > 4 that there is some value of k such that for all greater values of k the values 
of R/B^ :_1 increase. The order k must surely be greater than this value before it can be claimed that the 
asymptotic regions has been achieved. If we use that maximum ratios as obtained from Table IT51 we estimate 
that for D = 4 we have rj r h ~ 0.12 and for D — 5 we have rj r h ~ 0.052. In D = 3 if we assume that the data 
of Table El extrapolates to a constant as k — > 00 then the radius of convergence of the Ree-Hoover ring is 
rjrh = 0.25. At most the ratios are bounded below by 0.91 which leads to r\ r h = 0.27. All of these estimated 
radii of convergence are substantially less than the freezing densities 77/ . 

In order for the radius of convergence of the virial series to be larger than these estimates obtained from 
the Ree-Hoover ring there must be cancellations between diagrams in the class Bk[k,i]. Such cancellations 
can occur because, for example, the diagrams R^ffc] and Rs^fc] which have magnitudes comparable to R[fc] 
have signs opposite to R[k]. Furthermore from Table ^ we see that Bk is indeed less in magnitude than 
R[fe] for k = 6, 7,8. What cannot be inferred from the existing data is whether or not when k is sufficiently 
large that the diagrams are in their asymptotic region that the cancellation is severe enough to reduce the 
exponential rate of growth of the ratios Bk+i/Bk- Unless the cancellation becomes sufficiently great for the 
radius of convergence to be larger than the freezing density then the leading singularity cannot be on the 
real 77 axis and there must be oscillation in the signs of the virial coefficients. 



6 Approximate equations of state for hard spheres 

For over 40 years the eight virial coefficients of Table^have been used to inspire many approximate equations 
of state for the low density phase of hard spheres. These approximates may be grouped by the location of 
their leading singularity into the following three classes: 

1: High order poles at 77 = 1 

Examples of these are the proposals of Thiclc 19 

1 + 277 + 3r? 2 

Pv/k B T = — - — — , (30) 

Reiss, Frisch, and Lebowitz (SU|, and Wertheim [3U] 

Pv/k B T= \\ ■' 1 . cm 

Guggenheim [2~2"| 

Pv/k I} 1 (32i 

and the proposal of Carnahan and Starling [2H 

w = : 1 ( ; :; )3 - ■ (33) 

We do note that Torquato |32l I33| proposes an equation of state which agrees with Eq. [231 for 77 < 77/ but 
which is of a different form for 77 > 77/. 

2: Simple poles at |261 125| or near the packing fraction r/ cp = 0.74048 - of closest packed hard 
spheres P |2"1 IT51 EH] 

In particular the Pade analysis made in 13 has simple poles at 

77 = (1.22 ± 0.09i)rj cp (34) 

which leads to sign oscillations beginning with -B45. 



(1 


-v) 2 ' 


1 + 


7? + 77 2 


(1 


-77)3 ' 




1 


(1 


-77)4' 




4- rj 2 — 77' 


(1 


-T7) 3 
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3: A fractional power law divergence at or near the "random close packed" density i] rcp = 0.64 
as denned by |61M62ll5lfll64) 

These approximates are obtained from a D-log Pade analysis and are (generalizations) of the form 

Pv/k B T = A(rj - rjrcp)- 3 (35) 

As an example s is estimated as 1 in [37] as 0.678 in [2B] and 0.76 in [31]]. In [5S] other values of r] rcp are 
chosen and the values of s lie in the range 0.6 < s < 0.9 depending on the approximation used. 

All these approximate equations of state share the feature that their leading singularity is at a value of 
rj which is greater than the freezing density 77/ = 0.49 and indeed is even greater than the solid end of the 
phase transition r/ s = 0.542. This fact has lead to the assumption that the virial expansion is analytic at the 
freezing transition r\f and this qualitative feature is incorporated into most of the phenomenological theories 
used to describe freezing J3j)] [3J5] [3H] . It is therefore of great importance that the estimates of the radius 
of convergence made above on the basis of loose packed dominance and the assumption of no cancellation 
estimated the radius of convergence at no more than 0.27. 

The estimate of the radius of convergence relies on values of k beyond the first eight virial coefficients 
used to obtain the approximate equations of state. Therefore it is fair to say that none of the approximates 
incorporates the true large k behavior of the virial coefficients. Thus even if there is cancellation for large k for 
the set of close packed diagrams none of the approximates is based on computations which can observe these 
cancellations and hence no known approximate equation of state can be considered reliable. We therefore 
conclude that at present there exists no evidence to support the claim that the virial expansion has a radius 
of convergence greater than the freezing density 77/. 
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A Individual diagram contributions for B4, B5, and Bq. 

In this appendix we tabulate the contributions of the individual Ree-Hoover diagrams to the virial coefficients 
for B4, B 5 , and Bq. 



Table 12: Individual diagram contributions to B^. 



D 


Ba\Q 11 IB'i 


B<f4 11 BA 




3 


31 673(2) 


—0 029781 f8l 


2869495 ■ ■ • 


4 


0.1888655 ■ ■ ■ 


—0.0370195 ■ ■ ■ 


0.1518460 ■ ■ • 


5 


0.115211(3) 


-0.039233(3) 


0.075978(4) 


6 


0.0714700 ■ ■ ■ 


—0.0381069 ■ • ■ 


0.03336314 ■ • ■ 


7 


0.044927(2) 


-0.035055(3) 


0.009873(3) 


7.7 


0.032669(2) 


-0.032331(3) 


0.000338(3) 


7.8 


0.031227(2) 


-0.031920(3) 


-0.000693(3) 


8 


0.0285344 ■ ■ • 


-0.0310921 ■ • ■ 


-0.0025576 ■ ■ • 


9 


0.018286(1) 


-0.026861(3) 


-0.008575(3) 


10 


0.0117986 ■ ■ • 


-0.0227611 ■ ■ ■ 


-0.0109624 ■ ■ • 


11 


0.0076638(8) 


-0.018997(3) 


-0.011333(3) 


12 


0.0050018 ■ ■ • 


-0.0156721 ■ • ■ 


-0.0109624- ■ ■ 


13 


0.0032819(5) 


-0.012805(2) 


-0.009523(2) 


14 


0.0021615(4) 


-0.010381(2) 


-0.008220(2) 


15 


0.0014288(3) 


-0.008362(2) 


-0.006933(2) 


20 


0.00018830(6) 


-0.0026504(7) 


-0.0024621(7) 


25 


0.00002615(1) 


-0.0007841(3) 


-0.0007580(3) 


30 


3.763(3) x 10~ 6 


-0.0002233(1) 


-0.0002196(1) 


35 


5.560(7) x 10~ 7 


-0.00006217(3) 


-0.00006162(3) 


40 


8.38(1) x 10~ 8 


-0.00001705(1) 


-0.00001697(1) 


45 


1.284(3) x 10~ 8 


-4.631(4) x 10" 6 


-4.618(4) x 10" 6 


50 


1.992(5) x 10~ 9 


-1.249(1) x 10" 6 


-1.247(1) x 10" 6 



Table 13: Individual diagram contributions to B$. Values for discs and spheres taken from [§]. The contri- 
butions from the Ree-Hoover ring diagram B$[5, 2] are underlined. 





discs 


spheres 


D = 4 


D = 5 


B 5 /B| 


0.3336 


0.1103 


0.03565(5) 


0.01297(1) 


B 5 [0,l]/B« 


0.3618 


0.1422 


0.059015(9) 


0.025442(1) 


B 5 [4,l]/B* 


-0.0266 


-0.0314 


-0.02650(2) 


-0.019184(5) 


B 5 [5,l]/B| 
B 5 [5,2]/Bf 
B 5 [5,3]/B 9 4 


-0.0102 
0.0086 



-0.0165 
0.0162 
-0.0002 


-0.01762(4) 
0.02131(2) 
-0.0005498(5) 


-0.015511(4) 
0.022980(7) 
-0.0007622(3) 




D = 6 


D = 7 


D = 8 


D = 50 


B 5 /B* 


0.007528(8) 


0.007071(7) 


0.007429(6) 


2.17(1) x 10~ 8 


B S [0,1]/B| 


0.0112852(7) 


0.0051189(4) 


0.0023640(3) 


1.67(8) x 10- Ib 


B 5 [i,l\/B- 


-0.012899(4) 


-0.008296(3) 


-0.005185(2) 


-6.1(1) x 10 _la 


B s [5,l]/B| 
B 5 [5,2]/B| 
B 5 [5,3]/B 2 4 


-0.012351(3) 
0.022332(6) 
-0.0008395(4) 


-0.009220(3) 
0.020277(6) 
-0.0008090(4) 


-0.006588(3) 
0.017522(6) 
-0.0007149(4) 


-1.38(3) x 10" 11 
2.17(1) x 10~ 8 
-4.84(7) x 10 -11 
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Table 14: Individual diagram contributions to B§. Values for discs and spheres taken from The contri- 
butions from the Ree-Hoover ring diagram ^[6,3] are underlined. For D — 50 when the value x.xx x 10 m 
has m < —16 we write ~ 10 m . 





discs 


spheres 


D = 5 


D = 6 


Be/B 2 


0.1994 


0.0386 


0.00102(8) 


-0.00176(2) 


fl 6 [0,l]/B§ 


0.2292 


0.0588 


0.0048248(9) 


0.0014771(1) 


B 6 [4, l]/B* 


-0.0273 


-0.0212 


-0.00569(2) 


-0.002600(3) 


B 6 [5,l]/B| 


-0.0191 


-0.0187 


-0.00719(1) 


-0.003800(4) 


B 6 [5,2]/B| 


0.0090 


0.0099 


0.00498(2) 


0.003038(4) 


B 6 [5 : 3]/B = 





-0.0002 


-0.000244(1) 


-0.0001759(3) 


B 6 [6,l]/-B2 


0.0088 


0.0132 


0.01029(4) 


0.00735(1) 


B 6 [6,2]/B 2 B 


0.0077 


0.0121 


0.01064(4) 


0.008204(8) 


B 6 [6,3]/B| 


-0.0051 


-0.0109 


-0.01702(5) 


-0.01693(2) 


B 6 [6,4]/B| 


-0.0019 


-0.0027 


-0.001520(4) 


-0.0009332(8) 


B 6 [6,5]/B= 


-0.0010 


-0.0022 


-0.001559(9) 


-0.001001(2) 


B 6 [6,6]/B^ 


-0.0009 


-0.0011 


-0.000553(1) 


-0.0003335(3) 


B 6 [6,7]/B| 


-0.0005 


-0.0007 


-0.000437(1) 


-0.0002948(3) 


B 6 [6,8]/Bf 


0.0004 


0.0008 


0.000712(2) 


0.0004958(8) 


B 6 [6,9]/b5 


0.0001 


0.0012 


0.00256(2) 


0.002356(7) 


B 6 [6, 10] /B° 


0.0000 


0.0002 


0.0003086(8) 


0.0002481(2) 


B 6 [6, 11] /B| 


-0.0000 


-0.0003 


-0.0002596(7) 


-0.0001587(1) 


B 6 [6, 12] IB\ 





-0.0002 


-0.000472(2) 


-0.0003931(6) 


B 6 [6, 13] /B| 





0.0002 


0.000318(4) 


0.000285(1) 


B 6 [6, 14] /B| 





-0.0000 


-0.0000991(5) 


-0.0001007(2) 


B 6 [6, 15] IB% 





0.0000 


0.000046(2) 


0.0000486(9) 


B 6 [6, 16] /Bf 





0.0004 


0.00138(1) 


0.001463(2) 


B 6 [6, 17]/B.= 





-0? 


-2.00(1) x 10" 6 


-3.176(6) x 10" 6 


B 6 [6, 18] jB% 








2.59(1) x 10" 7 


3.696(8) x 10" 7 





D = 7 


D = 8 


D - 50 


B 6 /Bl 


-0.00352(2) 


-0.00451(1) 


-7.6(2) x 


10" 


10 


B 6 [0,l]/B b 2 


0.00046725(4) 


0.00015174(2) 


~ 1Q-'^ 


B 6 [4, 1]/B 2 S 


-0.001145(1) 


-0.0004946(6) 


~ 10 


B 6 [5,l]/B° 


-0.001902(2) 


-0.0009209(9) 


~ io" ls 






B 6 [5,2]/B| 


0.001752(2) 


0.0009748(9) 


~ io" 16 






B 6 [5,3]/B| 


-0.0001124(2) 


-0.0000665(1) 


nn -18 
~ 10 






Be [6, 1]/B 2 5 


0.004866(6) 


0.003063(3) 


6.7(7) x 


10" 


lb 


B 6 [6,2]/B| 


0.005887(5) 


0.004019(3) 


1.6(1) x 


10" 


13 


B 6 [6,3]/B| 


-0.01540(2) 


-0.01318(2) 


-7.6(2) x 


10" 


10 


B 6 [6,4]/B| 


-0.0005349(5) 


-0.0002926(3) 


~ 10 






B 6 [6,5]/B! 


-0.000589(1) 


-0.0003277(6) 


~ 10 






B 6 [6,6]/B= 


-0.0001883(1) 


-0.00010175(8) 


~ io" 18 






B 6 [6,7]/B 2 5 


-0.0001858(2) 


-0.0001114(2) 


~ io" 17 






B 6 [6,8]/B| 


0.0003158(4) 


0.0001895(2) 


~ io" 17 




14 


B e [6,9]/B% 


0.001912(4) 


0.001428(3) 


5.3(5) x 


10" 


B 6 [6, 10] /Bf 


0.0001778(1) 


0.0001180(1) 


~ IO" 16 






B 6 [6, ll]/B.f 


-0.00008782(7) 


-0.00004567(4) 


~ io" 20 






B 6 [6, 12] IB% 


-0.0002848(4) 


-0.0001889(2) 


~ io" 17 






B 6 [6, 13] /B| 


0.0002231(6) 


0.0001595(4) 


~ io- 16 






B 6 [6, 14] /B| 


-0.0000868(2) 


-0.0000671(1) 


in— 18 
~ 10 






B 6 [6, 15] IB% 


0.0000447(4) 


0.0000371(3) 


2.2(3) x 


10" 


15 


B 6 [6, 16] /Bf 


0.001358(2) 


0.001147(2) 


1.34(8) x 10 


-12 


B 6 [6, 17] /Bf 


-3.774(8) x 10" 6 


-3.758(9) x 10" 6 


-1.5(1) x 


10" 


15 


B 6 [6, 18] IB% 


3.944(9) x 10" 7 


3.524(9) x 10" 7 


nn -is 
~ 10 
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B Numerical values of selected diagrams to high order. 

In this appendix we tabulate the results of Monte-Carlo evaluations of selected diagrams Bk[m,i] to orders 
up to k = 17. 



Table 15: 0/5 2 fe ~ 1 = B k [0, 



k 


D = 2 


D = 3 


D = 4 


D = 5 


D = 6 


3 


0.7821(1) 


0.6248(2) 


0.5063(2) 


0.4143(2) 


0.3410(2) 


4 


0.5488(4) 


0.3166(3) 


0.1888(2) 


0.1153(2) 


0.0713(2) 


5 


0.3620(3) 


0.1420(2) 


0.0591(2) 


0.02522(8) 


0.01121(7) 


6 


0.2292(3) 


0.0593(2) 


0.01648(6) 


0.00487(6) 


0.00148(2) 


7 


0.1412(3) 


0.0233(2) 


0.00424(6) 


0.00076(1) 


0.000170(3) 


8 


0.0844(4) 


0.0087(2) 


0.00101(2) 


0.000129(3) 


1.81(5) x 10~ 5 


9 


0.0505(4) 


0.00315(6) 


0.000226(5) 


1.78(7) x 10~ 5 


1.3(2) x 10~ 6 


10 


0.0293(4) 


0.00111(2) 


5.2(2) x 10~ 5 


2.5(4) x 10~ 6 




11 


0.0170(3) 


0.000380(8) 


1.0(1) x 10~ 5 






12 


0.0097(2) 


0.000128(3) 


2.7(7) x 10~ 6 






13 


0.0053(1) 


5.2(4) x 10~ 5 








14 


0.00304(6) 


1.7(3) x 10~ 5 








15 


0.00179(4) 










16 


0.00098(2) 










17 


0.00055(1) 











Table 16: j [/B^ 1 = B fc [4, l]/B% 



k 


D = 2 


D = 3 


D = 4 


D = 5 


D = 6 


4 


-0.01644(5) 


-0.02981(9) 


-0.0370(1) 


-0.0391(1) 


-0.0382(1) 


5 


-0.0264(3) 


-0.0316(2) 


-0.0270(2) 


-0.0189(2) 


-0.0130(2) 


6 


-0.0285(5) 


-0.0219(4) 


-0.0117(2) 


-0.0059(1) 


-0.00256(5) 


7 


-0.0239(5) 


-0.0114(2) 


-0.00403(8) 


-0.00130(3) 


-0.00039(1) 


8 


-0.0183(4) 


-0.0056(1) 


-0.00120(5) 






9 


-0.0123(3) 


-0.0025(1) 








10 


-0.0086(4) 










11 


-0.0056(2) 










12 


-0.0038(2) 










13 


-0.0023(3) 











Table 17: /B^ 1 = B k [b, 1]/B 2 fe 



k 


D = 2 


D = 3 


D = 4 


D = 5 


D = 6 


5 


-0.01016(6) 


-0.0164(1) 


-0.01741(9) 


-0.01550(9) 


-0.0123(1) 


6 


-0.0188(3) 


-0.0189(2) 


-0.0124(2) 


-0.0069(1) 


-0.00386(8) 


7 


-0.0227(3) 


-0.0130(3) 


-0.0053(1) 


-0.00200(4) 


-0.00070(1) 


8 


-0.0205(4) 


-0.0072(1) 


-0.00199(5) 






9 


-0.0163(3) 


-0.0036(1) 








10 


-0.0115(5) 










11 


-0.0083(2) 










12 


-0.0056(3) 










13 


-0.0043(5) 
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Tabic 18: R/-B2 The underline marks the approximate location of the minimum value. 



k 


D = 2 


D = 3 


D = 4 


3 


0.7824(2) 


0.6248(2) 


0.5063(2) 


4 


— u.uiooy(y j 


n none 1 1 \ 

— u.uzyo(i ) 


— U.Uo ( 1(1) 


5 


0.00860(6) 


0.01623(9) 


0.0214(2) 


g 


—0 00526(81 


—0 0109(11 


—0 0150(21 


7 


00335(6") 


0078(21 


0.0124(2) 


8 


—0 00234(51 


—0 0064(11 


—0 0106(21 


9 


001 77(4") 


0053(11 


0098(21 


10 


—0 001 25(31 


—0 00452(91 


—0 0091 (21 


11 


0.00095(2) 


0.00392(8) 


0.0089(2) 


12 


-0.00074(1) 


-0.00333(7) 


-0.0083(2) 


13 


0.00055(1) 


0.00313(8) 


0.0086(2) 


14 


-0.00041(1) 


-0.0027(1) 


-0.0086(2) 


15 


0.00033(2) 


0.0026(1) 


0.0087(3) 


16 


-0.00023(2) 






17 


0.00021(3) 






k 


D = 5 


D = 6 


D = 7 


3 


0.4139(2) 


0.3409(2) 


0.2822(2) 


4 


-0.03925(5) 


-0.03815(8) 


-0.0351(1) 


5 


0.0230(1) 


0.02210(6) 


0.02025(9) 


6 


-0.01689(9) 


-0.0173(1) 


-0.0153(2) 


7 


0.0142(3) 


0.0144(2) 


0.0135(2) 


8 


-0.0129(3) 


-0.0141(3) 


-0.0132(2) 


9 


0.0126(2) 


0.0138(3) 


0.0132(3) 


10 


-0.0126(3) 


-0.0150(3) 


-0.0143(3) 


11 


0.0128(3) 


0.0162(3) 




12 


-0.0134(3) 


-0.0172(3) 




13 


0.0142(3) 


0.0201(4) 




14 


-0.0166(3) 


-0.0238(5) 




15 


0.0183(4) 


0.0281(6) 





Table 19: R (-^-) /B* 1 = ^4,2/^3 1 ■ The underline marks the approximate location of the minimum value. 



k 


D = 2 


D = 3 


D = 4 


D = 5 


D = 6 


5 


-0.01017(4) 


-0.01654(7) 


-0.01748(6) 


-0.01551(4) 


-0.01235(3) 


6 


0.00756(6) 


0.01191(8) 


0.01257(6) 


0.01067(6) 


0.00816(4) 


7 


-0.00607(8) 


-0.01003(7) 


-0.0106(1) 


-0.00936(9) 


-0.00737(7) 


8 


0.0051(1) 


0.0088(2) 


0.0100(1) 


0.00902(8) 


0.00722(7) 


9 


-0.00427(8) 


-0.0083(1) 


-0.0100(1) 


-0.0095(2) 


-0.00762(8) 


10 


0.00381(8) 


0.0082(2) 


0.0106(2) 


0.0106(2) 




11 


-0.00309(3) 


-0.0075(1) 


-0.0109(2) 


-0.0116(2) 




12 


0.00259(3) 


0.00728(7) 


0.0115(2) 


0.0136(3) 




13 


-0.00220(3) 


-0.0071(3) 


-0.0125(3) 


-0.0156(3) 




14 


0.00188(4) 


0.0068(4) 


0.0140(5) 


0.0172(4) 




15 


-0.00151(5) 


-0.0060(5) 


-0.0145(7) 


-0.0227(7) 




16 


0.00131(6) 










17 


-0.00096(8) 
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Table 20: R (-O - ) 7-^2 1 = ^4,1/^2 1 ■ The underline marks the approximate location of the minimum value. 



k 


D = 2 


D = 3 


D = 4 


D = 5 


5 




-0.000242(2) 


-0.000550(2) 


-0.000765(4) 


6 




0.000435(3) 


0.000984(7) 


0.00137(1) 


7 


-1.03(5) x 10~ s 


-0.000342(5) 


-0.00078(1) 


-0.00112(1) 


8 


5.0(2) x 10~ 7 


0.000294(6) 


0.00076(1) 


0.00107(1) 


9 


-3.5(2) x 10~ 6 


-0.000290(6) 


-0.00073(1) 


-0.00108(2) 


10 


8.6(4) x 10~ 6 


0.000278(6) 


0.00074(1) 


0.00121(2) 


11 


-1.53(8) x 10" 5 


-0.000281(6) 


-0.00082(2) 


-0.00135(3) 


12 


1.83(9) x 10~ 5 


0.000273(8) 


0.00085(2) 


0.00154(3) 


13 


-1.90(9) x 10~ 5 


-0.00029(1) 


-0.00090(3) 


-0.00178(4) 


14 


2.1(1) x 10~ 5 


0.00028(2) 


0.00104(4) 


0.00205(5) 


15 


-1.79(9) x 10" 5 


-0.00026(2) 


-0.00102(6) 


-0.00244(9) 


16 


1.6(1) x 10~ 5 








17 


-1.3(1) x 10~ 5 








18 


1.3(2) x 10" 5 









Table 21: R {^/B^ 1 = R^/B^ 1 



k 


D = 2 


D = 3 


D = 4 


6 


0.00848(7) 


0.01328(9) 


0.01286(5) 


7 


-0.0063(1) 


-0.0094(1) 


-0.00947(9) 


8 


0.00469(9) 


0.0078(1) 


0.0081(1) 


9 


-0.00363(7) 


-0.0065(1) 


-0.00748(9) 


10 


0.00290(6) 


0.0061(1) 


0.0076(1) 


11 


-0.00230(8) 


-0.0056(1) 


-0.0072(1) 


12 


0.0020(1) 


0.0051(1) 


0.0075(2) 


13 


-0.0017(2) 


-0.0045(2) 


-0.0080(2) 


14 


0.0013(2) 


0.0048(3) 


0.0077(3) 


15 


-0.00097(24) 


-0.0043(4) 


-0.0087(5) 



Table 22: Pinwheel/B^ 



k 


D = 2 


D = 3 


D = 4 


5 


-0.0266(1) 


-0.0315(1) 


-0.02647(6) 


6 


0.00906(9) 


0.01000(8) 


0.00754(6) 


7 


-0.00194(4) 


-0.00230(3) 


-0.00178(3) 


8 


0.000187(4) 


0.000302(6) 


0.000261(5) 


9 


-3.7(1) x 10~ 6 


-0.0000157(5) 


-0.000022(1) 


10 


2(1) x 10- 8 


5(1) x 10~ 7 


8(4) x 10- 7 
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C Graphs of the ratio of selected diagrams to high order. 

In this appendix we graph the ratios -Bfc[4, 1]/Bk[0, 1], Bk[5, 1]/Bk[0, 1], R4,2/R, R4,i/R, and Rs^/R to orders 
up to k = 17. 



1.5 



D=4 



B k [4,l]/B k [0,l] 



D=3 



0.5 - 



D=2 



10 

k 



15 



Figure 1: Absolute value of j/0 = Bfe[4, 1]/Sfe[0, 1] in dimensions 2 (triangles), 3 (filled circles), and 4 
(crosses). 



2.5 



2 - 



D=4 



B fc [5,l]/B fc [0,l] 



1.5 - 



0.5 - 



D = 3 



D = 2; 



10 

k 



15 



Figure 2: Absolute value of /0 = 2?fc[5, 1]/Bfc[0, 1] in dimensions 2 (triangles), 3 (filled circles), and 4 
(crosses). 
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D = 2 



R4 ;2 /R 



i i 



i . D = 3 



D = 4 
□ D = 5 



10 



15 



Figure 3: Absolute value of R(-^>-)/i? = R4,2/R in dimensions 2 (triangles), 3 (filled circles), 4 (crosses), and 
5 (squares). 



R 4 ,i/R 



0.15 



0.1 



0.05 - 



1 



s * # 



D = 5 
D = 4 
D = 3 



D=2- 



10 



15 



Figure 4: Absolute value of R(-Q-)/R = R4 j i/R in dimensions 2 (triangles), 3 (filled circles), 4 (crosses), and 
5 (squares). 



21 



: } 






, D = 2 - 


















- * 


i 


d=3 ; 


* 














• • ■ * 




i D = 4 1 



k 

Figure 5: Absolute value of R(-^-)/R = Rs^/R in dimensions 2 (triangles), 3 (filled circles), and 4 (crosses). 
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